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The exactly solvable four-vertex model on a square grid with the different bound- 
ary conditions is considered. The application of the Algebraic Bethe Ansatz method 
allows to calculate the partition function of the model. For the fixed boundary con- 
ditions the connection of the scalar product of the state vectors with the generating 
function of the column and row strict boxed plane partitions is established. Tiling 
model on a periodic grid is discussed. 

I. INTRODUCTION 

The connection of integrable models with the enumerative combinatorics is widely dis- 
cussed. The six- vertex model plays an important role in these considerations. The model was 
widely studied both for periodic boundary conditions and for the fixed ones [l[ 2], \a, O, Q]. 
The six-vertex model for the so-called fixed boundary conditions is related to the enumera- 
tion of domino tilings of Aztec diamonds and of alternating sign matrices 0, Isl, EJ- 

the six- vertex model in which two vertices 



11| where the connection of this model 



The four-vertex model is a particular case o: 
are frozen out. This model was considered in lo|, 
with a random tiling models on a semi-infinite strip geometry was discussed. Tiling models 
are classical discrete statistical models in which the tiles of the different geometric shapes 
are packed together so that they cover space without holes or overlaps. They are called 
random only to emphasize the difference with perfect quasi-periodic tilings. 

In this paper we shall solve the model by the direct application of the Quantum Inverse 
Scattering Method (QISM) Q, Q. This approach will allow us to study the model on a 
finite lattice with the different boundary conditions and to calculate partition function of 
both vertex and corresponding tiling models. In particular, we shall show that the partition 
function of the four-vertex model with the fixed boundary conditions is related to the gen- 
erating function of the boxed plane partitions. The boxed plane partitions are equivalent to 
the three dimensional Young diagrams placed into a box of a finite size and to lozenge tiling 



of an semiregular hexagon 

model is equivalent to the lozenge tiling of the torus 
The connection of the boson type integrable mode 



For the periodic boundary conditions the four-vertex 
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and with the theory of plane partitions 



s with the theory of symmetric 



161 ] was considered in the papers 



The paper is organized as follows. We give the definition of the model and discuss the 
QISM approach to its solution in Sec. II. In Sec. Ill the model is studied on the lattice with 
the fixed boundary conditions and its connection with the plane partitions is established. 
The partition function of the model on the periodic grid is calculated in Sec. IV. 



II. FOUR- VERTEX MODEL 



A special case of the six vertex model - a four vertex model on a square grid is described 
by four different arrows arrangements pointing in and out of each vertex. A statistical 
weight corresponds to each type of the vertices and there are three vertex weights u a , u b 
and uj c (FIG. 1). There is an alternate description of the vertices in terms of lines flowing 
through the vertices. Since a lattice edge can exist in two states, line or no line, there exists 
a one-to-one correspondence between the arrow configuration on the lattice and the graphs 
of lines on the lattice - nests of lattice paths. 
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FIG. 1: The four allowed types of vertices in terms of arrows and lines. 

For the homogeneous model the vertex weights do not depend on the position of the 
vertex and the partition function is equal to 

Z(u a , U b , U) c ) = J2 WaW l b U l °, (1) 

where the summation is extended over all allowed configurations of the arrows on a lattice. 
These configurations depend on the imposed boundary conditions which are specified by 
the direction of the arrows on the boundary of the grid. The number of vertices (a), (b), (c) 



in each configuration is l a ,l b ,l c . For the inhomogeneous model the vertex weights are site 
dependent. 

To apply the Quantum Inverse Scattering Method we use the spin description of the 
model. With each vertical bond of the grid and horizontal bond one associates the space 
C 2 , with spin up and spin down states forming a natural basis in this space. The spin up 
state on the vertical bond corresponds to the line pointing up, while the spin down state 

(A 

to the line pointing down. The spin up state on the z-th horizontal bond =1 *~ )i 

corresponds to the horizontal line pointing to the left, spin down state to the line pointing 
to the right I 1=1 The total space of the vertical lines is V = (C 2 ) 0JV and we 

shall call it an auxiliary space. The total space of the horizontal lines is Ti = (C 2 )® N , and 
we shall call it a quantum space. With each vertex of the lattice one associates an operator 
acting in the full space V <S> 7i. This operator is called L-operator and it acts nontrivially 
only in a single vertical space C 2 and in single horizontal space C 2 , while in all other spaces 
it acts as the unity operator. 

The L-operator of the four vertex model is equal to 

L(n\u) = --(1 + a z )e n + — (1 - a z )e n + a + a~ + <r~a+, (2) 

where a z,± - are the Pauli matrices, the projector e = ^(cr z + l), and the matrix with subindex 
n acts nontrivially only in the n-th space: s n — I ® ... ® I ® s ® I ® ... ® I . 
One can represent the matrix elements of the introduced L-operator 

/ 



L(n\u) 



(3) 



L u (n\u) L l2 {n\u) \ _ / -ue n a n 
\L 2 i(n\u) L 2 2{n\u) J I cr+ u^e 

as a dot with the attached arrows (see FIG. 2). The matrix element Lu(n\u) corresponds to a 
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FIG. 2: Vertex representation of the matrix elements of the L-operator. 



vertex (1) (FIG. 2), where a dot stands for the operator — ue n , this operator acts on the local 



spin state and the only non-zero matrix element of this operator is n (<— \ — ue n \ <—) n what 
gives the vertex (b) (FIG. 1) with a weight — u. The matrix element L 2 2{n\u) corresponds 
to a vertex (2) (FIG. 2), where a dot is the operator u~ x e n and the non-zero matrix element 



\u 



-l. 



-) n is the vertex (a) (FIG. 1) with a weight u 1 . The matrix elements Li2(n\u), 



and Li2\(n\u) correspond to the vertices (3) and (4) (FIG. 2) respectively, and the nonzero 
matrix elements n {— > \a~\ <—) n and n {<— |cr+| — ►)„ are the vertices (c) (FIG. 1) with a weight 
1. 

The L-operator §2§ satisfies the intertwining relation 

R(u, v) (L(n\u) ® L(n\v)) = (L(n\v) ® L{n\u)) R(u, v), 



(4) 



in which R(u, v) is the 4x4 matrix 

/ 



R(u, v) 
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(6) 



The vertical monodromy matrix is the product of L-operators 

/ 

T{u) = L(M\u)L(M - l\u)...L(0\u) = 



(7) 



A{u) B(u 
\C(u) D(u 

The matrix entries of the monodromy matrix (J7D are expressed as sums over all possible 
configurations of arrows with different boundary conditions on a one-dimensional lattice 
with M + l sites (FIG. 3). Namely, operator B{u) = Y* M ,... )kl=1 L lkM (M\u)L kMkM _ 1 (M - 
l\u) . . . Ljt 1 2(0|tt) corresponds to the boundary conditions when arrows on the top and bot- 
tom of the lattice are pointing outwards (configuration (B)). Operator C(u) corresponds to 
the boundary conditions when arrows on the top and bottom of the lattice are pointing in- 
ward (configuration (C)). Operators A(u) and D(u) correspond to the boundary conditions 
when arrows on the top and bottom of the lattice are pointing up and down respectively 
(configurations (A) and (D)). 

The commutation relations of the matrix elements of the monodromy matrix are given 
by the same R— matrix (JSj) 

R(u, v) (T(u) ® T(v)) = (T(v) ® T(u)) R(u, v). (8) 



(B)(C)(A)(D) 



FIG. 3: Graphic representation of the entries of the monodromy matrix 

The most important relations are 

C(u)B(v) = g(u, v) {A(u)D(v) - A(v)D(u)} , (9) 

A(u)B(v) = f{u,v)B(v)A(u) + g{v,u)B(u)A(v), 

D(u)B(v) = f{v,u)B(v)D(u) + g{u,v)B(u)D(v), 

[B(u),B(v)\ = [C(u),C(v)]=0. 

The transfer matrix r(u) is the trace of the monodromy matrix in the auxiliary space 

t{u) = trT(u) = A{u) + D(u). (10) 

The relation ([H]) means that [t(u),t(v)] = for arbitrary values of parameters u,v. 
The L-operator ([2]) satisfies the relation 

e c <L(n\u)e- c < = e"^ L(n\u)e* a \ (11) 

where C is an arbitrary parameter. From this equation and from the definition of the 
monodromy matrix (I7j) it follows that 

e^ sz T{u)e^ sz = e~^T{u)e^\ (12) 

where S z = \ Y^iLq of is the operator of z component of the total spin. The equalities 

S z B(u) = B{u) {S z - 1) , (13) 
S*C(u) = C(u) {S z + 1) , 

are the consequence of the equation f|T2l . It means that the operator B{u) decreases the total 
spin of the system, while C(u) increases it. From the equation fll2p follows the commutativity 

[t(u),S*] = 0. (14) 



The generating state in the space Tt is the state with all spins up 



=> = ®&,| <->* = ®£o I 1 I • (15) 



It is annihilated by the operator C(u) 

C(u)| <«=) = <>, (16) 
and it is an eigenvector of operators A(u) and D(u), 

A(u)\ <=> = a(«)| -<=); £>(u)| <=) = <*=), (17) 

with the eigenvalues 

a(u) = {iu) M+ \ 5{u) = {-) M+1 . (18) 



it 



The total spin of the generating state is \{M + 1): 



S*|<=) = i(M+l)|<=). (19) 

One of the main objects in the further consideration will be the vector generated by the 
multiple action of operators B(u) on the state | <=) 

N 

\y n { Ul ,u 2 ,...,u N )} =Y[B(ui)\ <«=). (20) 

i=l 

From (1131) it follows that the total spin of this vector is 

N i N 

S z n B( Ul )\ -<=) = - (M + 1 — 2N) JJ S( Wl )| <=). (21) 

i=l i=l 



The state conjugated to (1201) is 

v 

i=l 

It is easy to verify that (<= \B{u) = 0. 

In the coordinate representation the state vector takes the form 

N 

\^n{ui,u 2 , ..,%)) = n B i u k)\ <=) = J2x fl (u 1 ,u 2 , ...jU^lmi, ...,m n ). (23) 

fc=i 

Here we denote the state with iV spins down in the sites mi, ...,mjy and with M + 1 — iV 
spins up in the other sites by |mi, mjv). The wave function Xn{ u u u 2i ■■•i u n) satisfies the 



exclusion condition. It is non equal to zero only if numbers ji = (m 1 , m 2 , rajy) form a 
strict partition M > m\ > m 2 > ... > m N > and mi > m i+1 + 2. By the direct calculations 
it may be shown that the wave function is of determinantal form: 

det ( u ~ 2 ( m *+ k ~ N )) 

U2, -..,«*) = (u^.U^' 2 ^ - ^ ^) 

lll<k<j<N\ u k ~~ U j ) 

In the terms of the symmetric Schur function 

N-k+\ k 



S\(ui,u 2 , ...,u N ) 



det ( Uj 



Yil<k<j<N( u k — Uj) 

it has the form 

X,( Ul ,u 2 ,...,u N ) = (-1)5>* ( Ul ... UN ) M - 2 ( N -V s^_ s (uf,uf,...,u- N 2 ) , (25) 

where 5 is a partition with the elements 5j = 2(N — j), j = 1, N, and the elements of 
the partition A = fi — 5 are Xj = rrij — 2(N — j) satisfy the condition M — 2(N — 1) > Ai > 
A 2 > ... > X N > 0. 



III. FIXED BOUNDARY CONDITIONS 

Let us consider now a model on a 2N x (M + 1) square lattice with a following boundary 
conditions: all arrows on the left and right boundaries are pointing to the left while the 
arrows on the top and bottom of the first iV columns (counting from the left) are pointing 
inwards and the arrows on the top and bottom of the last iV ones are pointing outwards. 
We shall call this condition - the fixed boundary condition. 

To enumerate all possible configuration of the vertices it is more convenient to use the 
description of the model in terms of flowing lines and to represent the allowed configurations 
as the nests of lattice paths. The path is running from one of the iV down left vertices to 
the top iV right ones and always moves east or north. The paths cannot touch each other 
and arbitrary number of consequent steps are allowed in vertical direction while only one 
step at a time is allowed in the horizontal one. If the first columns have the numbers 
(-N, . . . , — 1) and the last ones (1, . . . , JV), if the lower row has the number and the upper 
one M, then the m-th path is running from the vertex (— N + m — 1; 0) to the vertex (m; M), 
1 < m < N. A typical nest of lattice paths is represented in (FIG. 4). The length of the path 
is N + M. The number of the (c) vertices in the admissible path is 2N since only one step is 



allowed in horizontal direction, so the number of (b) vertices is M — N + l. It means that the 
number of (c) and (b) vertices in the nest is equal to l c = 2N 2 and l h = l a = N(M — N + l) 
respectively. 
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FIG. 4: A typical configuration of admissible lattice paths with the fixed boundary conditions and 
the corresponding plane partition with the gradient lines. 



The partition function of the homogenous model with the fixed boundary conditions is 
equal to 

Z(u a ,u b ,u; c ) = (u a u b ) N ^ M - N+1 ^ 2 c N2 S(N,M), (26) 

where S(N, M) is the total number of the allowed nests of lattice paths. 

In the inhomogeneous model the vertex weights u a , Ub and uj c are site dependent. We 
shall consider a case when the weights depend on the number of the column only. If the first 
N columns have the numbers (-N, . . . , —1) and the last ones (1, . . . , N) then the partition 
function is equal to 



—N 



N 



k=-l j=l 

Consider the scalar product of the state vectors (T2"U]) and (T2"2]) 

W(ui,...,u N ;v u ...,v N ) = (<= \C(v 1 )...C(v N )B(u 1 )...B(u N )\ <=), 



(27) 



(28) 



where {u} and {v} are the sets of independent parameters. We can represent the matrix 
element ( !28l) as the two-dimensional square lattice with 2N x (M+ 1) sites. First iV vertical 
rows of the lattice are associated with the operators C(vj) and the last N vertical rows 
with operators B(uj). The horizontal rows of the lattice are associated with the local spin 
spaces, i-th row with the z-th space respectively. From the graphical representation of the 



operators B(u) and C(u) it follows that the matrix element (1281) is equal to the sum over 
all allowed configurations of vertices on a square lattice with the arrows on first N vertical 
rows pointing inwards, on the last iV ones pointing outwards; on the right and on the left 
boundaries all arrows are pointing to the left (FIG. 4): 

-N N 

W( Ul ,...,u N ; Vl ,...,v N ) = E II Ift-^W' ( 29 ) 

fc=-l j=l 



■iri n 

fe=-i j=i 



If we put 



{u Ja ) j = vZ), {w h )j = v-j] —N < j < -1, (30) 
(u a )j = uj 1 , (u h )j = Uj] 1 < j < N, 



(U c )j = 1 



then 



Z({u a },{u b },{u c }) = (-l) MN W( Ul ,...,u N ;v u ...,v N ). (31) 



The matrix element (1251) for the arbit rary values of the parameters Uj , Vj is evaluated by 



means of the commutation relations 



13J, 



171 ] and is equal to: 

W{u x ,...,u n] v u ...,v n ) = {-l) MN \\{^- 2 II -^\}detH, (32) 



j>k n - v] c> m n - u2 m 



where the entries of N x N matrix H are 



{/ x M-N+2 f s. -M+N-2) 
ft) "ft) 

Formulas ( 1321) . ( l33l) give the solution of the inhomogeneous model ( 1291) and represent its 
partition function in the determinantal form. 

The explicit answer for the partition function may be obtained if we choose the weight 
of the vertices 

Vj = q~*, uj = q 1 ^, (34) 

where q = e~ u , with v > is a Boltzmann weight. Under this parametrization the partition 
function ( 1271) will take a form 

Z(q) = £ qXjZ-i SPS-^+E^ ^C?^), (35) 



where the summation is over all admissible nests of lattice paths. The parametrization 
transforms (13"21) into 

W q (N,M) = (-1)™+^ |n (q** -?-^)" 2 JdetQ, (36) 
with the matrix elements 

k+j-1 k+j-1 

Qjk = k+j-1 k+j-1 > (37) 

q 2 -q 2 

and s = q M ~ N + 2 , The determinant of the matrix Q was calculated in |7| and is equal to 

!\ 1 j k 1 j ~ k 

[](g 2 -g 2 W II — ^Ei nsET-- ( 38 ) 
i>fc I i<j,k<N q 2 - q 2 



The substitution of this determinant into (!3"oT) gives 

W,(iV,M) = (_l)^V^+2-2iV) -Q 

l<j,k<N 



I — qM-N+2+j-k 
1 - qk+j-i 



N 2, 1 _ n M+3-j-k 
' i\NM -^ r {M+2-2N) TT 1 g 

~ A J 9 11 1 _ a k+j-l 

i<j,k<N 1 y 



and the partition function flHol) is given by 



w 2 , 1 _ n M+3-j-k 

z(s) = n .V,-, • (39) 

l<j,k<N " 

From the definition Z(l) = S(N,M) and for the partition function of the homogeneous 
model f[26"j) we have the following expression 

z( WttI ^, Wc )=( Wa ^^+ i ) W f 9 n M ~ J ~- 1 +3 - ^ 4 °) 

l<j,k<N J + K 1 

To connect the four-vertex model with the tiling model notice that each admissible con- 
figuration of lattice paths may be associated with an N x N array 7Ty. The m— th path 
may be thought of as the m— th column in this array with a matrix elements 7Tj TO equal to 
the number of the cells in the subsequent columns j of the lattice (starting from the right) 
under the m-th path. The array 

1 6 4 3^ 

TT = 5 3 1 (41) 

V 2 1 °J 

corresponds to the nest of paths in (FIG. 4). 



1 2 3 



FIG. 5: Three types of lozenges. 



An array 71"^- of non-negative integers that are nonincreasing as functions of both i and j 
(i,j = 1,2, ...) is called a plane partition 7r [14j], the integers 7Ty are the parts of the plane 
partition. Each plane partition has a three dimensional diagram which can be interpreted 
as stacks of unit cubes - the three-dimensional Young diagram, \n\ = being its volume. 
The height of a stack with coordinates is equal to the part of plane partition tt^-. If 
we have i < r, j < s and nij < t for all cubes of the plane partition, it is said that the plane 
partition is contained in a box with side lengths r, s, t. If tt^- > , i.e. if the parts of 

plane partition n are decaying along each column, then ir is called a column strict plane 
partition. We shall call the strict plane partition the partition n that is decaying along each 
column and each row (7r.y > 7Tj+ij and iTij > TCij+i). The element 7Tn of the strict plane 
partition 7r^ satisfies the condition Tin > 2r — 2 if all i,j < r. An arbitrary plane partition 
in a box r xr xt may be transferred into a strict plane partition in a box r x r x (t + 2r — 2) 
by adding to an array (71"^) the matrix 



7T 



2r - 2 2r - 3 • • • r - 1 
2r-3 2r-4 ••■ r-2 



r — 1 r — 2 







v ' - 1 u / 

which corresponds to a minimal strict plane partition. 

A plane partition in a r x s x t box is equivalent to a lozenge tiling of an (r, s, £)-semiregular 
hexagon. The term lozenge is referred to a unit rhombi with angles of | and ^ (FIG. 5). The 
lozenge tilings are simply the projections of three-dimensional diagrams with gradient lines. 
The three dimensional diagram which corresponds to the plane partition fj4TT) is represented 
in (FIG. 4). 

The partition function of the three dimensional Young diagrams placed into a box is equal 

to 

Z YD {q) = Y.^\ (42) 
where q is a Boltzmann weight and summation is over all plane partitions in a box. The 



volume of the three dimensional strict Young diagram in a box NxNxM which corresponds 
to the allowed arrows (paths) configuration is equal to [171 ] : 

, N 2 M ^ k A 



7T 



fc=-l " 3=1 

Substituting this expression into (1531) we obtain the partition function of the three dimen- 
sional strict Young diagrams 

.,->.. ]_ _ gM+3-j-k 



i<j,k<N 1 y 



(43) 



This partition function is the generating function of the strict plane partitions. The number 
of the correspondent lozenge tilings is Zyd{1) and is equal to the partition function of the 
homogeneous model ( 1261) with all weights equal to one: 

M - j - k + 3 



Z YD (1) = Z(l,l,l) 



n 



l<j,k<N 



j + k-1 



IV. PERIODIC BOUNDARY CONDITIONS 



Let us consider now the model on a lattice with the periodic boundary conditions. It 
means that on a square lattice with (M+l) horizontal rows and L vertical ones the boundary 
arrows on any horizontal row and on any vertical one are pointing in the same direction (see 
FIG. 6). Here we shall consider the case when (M + 1) and L are even. In the four- vertex 



M- 




FIG. 6: A typical configuration of arrows and lattice paths for the periodic boundary conditions. 



model the number of arrows pointing to the right or to the left is conserved in the subsequent 
columns of the grid. It means that the number of vertices c per column is also conserved 
and l c = 2n, where n is the number of arrows pointing to the right in a column. 



The partition function ([T]) of the homogeneous model on a periodic (M + 1) x L lattice 
with a fixed number of arrows pointing to the right and the vertex weights chosen to be 
ui a = e~ u , Ub = e u and u c = 1 is equal to 

ZM = Y,<? {h ~ la \ (44) 

where the summation is carried out over all allowed configurations of arrows on a lattice 
with the periodic boundary conditions or over all nests of the admissible lattice paths (FIG. 
6). 

Consider the trace of the transfer matrix ffTOl): 



trT L (u) = tr (A(u) + D{u)) L = ^(mj, m n \ (A(u) + D{u)) L \m u m n ). (45) 

Here we denote the state with n spins down in the sites mi, ...,m n and with M + 1 — n 
spins up in the other sites by |mi, ...,m n ), and the sum is taken over the complete set of 
states |mi,...,m n ) in the quantum space Tt. The number n is conserved because of the 
commutativity (1141) . Since t l (u) is a generating function of different combinations of A and 
D operators, the trace ( T45j) is the sum over all allowed configurations of vertices (a), (b), (c) 
with the weights S(u), a(u), 1 (|T8l) on a periodic (M + 1) x L grid 

trr L (u) = J2(*u) lb (in- 1 ) 1 " ■ (46) 

Comparing this expression with (jHJ) we obtain the following representation for the partition 
function 

Z n {v) = e-^ M+l ' 2n ^ L trr L (-ie u ), (47) 

where we put u = —ie u . 

To calculate this expression one has to solve the eigenvalue problem for the transfer 
matrix t{u) which means finding such a vector \§! n {y\, ...,v n )) that 

T(u)\V n (v 1 ,...,V n )) = Q n (u;Vi,...,V n )\^ n (v 1 ,...,V n )), (48) 

S z \V n (v u ...,v n )) = i(M + l-2n)|* n (vi, ...,!;„)). 

The vector \^ n (vi, ...,v n )) is an eigenvector of the transfer matrix if parameters Vj satisfy 
Bethe equations 



n 



1) V , k = l,...,n. (49) 

2 n 2 



v- 2 = n«j 



By putting v\ = e %Pk we can express solutions of the Bethe equations (jUJ) in the form 

2vr/ fc - P 

Pk = 1 <Pk < -^ (50) 

where /„ are integers or half-integers depending on n being odd or even, the total momentum 



P = J2]=iPj- These solutions were classified in [2lL l22|. 



The eigenvalue of the transfer matrix is equal to 

n -i 

9 n (n; {v}) = e <iC"+i-*0 {u M+1 V~ 2 + (-1)^-( M+1 - 2 ™)} J] — 2- (51) 

j= i u v j 

The partition function of the four vertex model with the periodic boundary conditions (|44|) 
is then given by the equation 

ZJU) = e -f (M+l-2n)X £ QL ( _ ie v. {v}) ^ (52) 

where summation is taken over all different sets of solutions of Bethe equations. The simplest 
case corresponds to n — 0. The partition function is Z Q (u) = {2cosh(M + l)v} L , and 
Zq{0) = 2 L . The other limiting case is when 2n = M + 1. In this case we have only two sets 
of solutions. The first is defined by the solutions (1501) of Bethe equations with Ik = k — 
k — 1, n and P = 0. The other one is given by = — ^r, k — 1, n, and P = — 7r. 
Applying the equality 

n 

IKa-be^) = a n -b n , 

k=l 

one can find that Z m+i (u) = 2. 

The vector conjugated to the Bethe eigenvector ( |48l) 

n 

\* n (p u ...,p n )) = l[B(^)\^), 

k=l 

where pu are the solutions (loUl) of Bethe equations (pfUl) . is equal to 

^ n ( Pl ,...,p n )\ = (^\f[B\e^). 

k=l 

The operators B(e ip ) and C(e ip ) are in the involution: B\e ip ) = (-l) M C(e ip ). 

The norm of the Bethe vectors is calculated with the help of ( |32i) . When the sets {u} = 
{v} the diagonal elements of the matrix H are equal to 



H H — M — n + 2. 



When {u} = {v} are the solutions of Bethe equations (1491) the nondiagonal entries of matrix 
H are units: Hjk = 1. It is not difficult to calculate the determinant 

det H = (M + 1)(M - n + l) ft_1 , 

and for the norm of any Bethe eigenvector we obtain: 



AT 2 ( Pl , ..., Pn ) = (* n (p!, ..., Pn )\V n ( Pl , ..., Pn )) = (M + 1)(M -n+ 1Y 1 - 1 n - 



e i Pj _ e i Pk 



where pj are the solutions (IqT?1) . The Bethe vectors (|48|) form a complete orthogonal set 22]. 

Consider the cylinder of circumference M + 1 and the length L. The probability that the 
lattice paths will enter the cylinder at the cites p, = (mi, rh n ); rhi > fh i+ i + 2 and leave it 
at the cites /i = (mi, m n ); m, > m i+ i + 2 is equal to: 

(mi, ...,m n |r L (-ze ly )|mi, ...,m n ) (53) 

= A/"~ 2 ^(mi, ...,m n |* n (pi, ...,p„))6^(-ie !/ ; {e ip })(^ n (pi, ...,p n )|mi, ...,m n ) 

{p} 

= M- 2 S,-s (e ip \e^, e*») (e~ ipi , e" ip2 , e"^) e^(-ie y ; {e*}), 
M 

the summation is taken over all solutions of Bethe equations, and representation ( |25l) has 
been used. 

We can associate the vertical and horizontal edges of the lattice carrying paths with 
lozenges. Lozenge (1) in (FIG. 5) corresponds to a vertical line of the path, while lozenge 
(2) to a horizontal one. Lattice edges without the paths correspond to lozenge (3). The 
number of lozenge tilings of the torus is given by (|52p with v — 0. 



V. CONCLUSION 

Though the considered four- vertex model is a particular case of the six-vertex model the 
determinantal representation of the partition function of the model with the fixed boundary 
conditions could not be obtained as a limit because of the absence of the correspondent 
answer for the six vertex case. The same is true for the representation of the mean value 
(153p in terms of the symmetric functions. 

The quantum Hamiltonian commuting with the transfer matrix of the four-vertex model 



is the Hamiltonian of the XX Z Heisenberg chain with the infinite anisotropy 



21 



22] . The 



approach described in this paper will allow to obtain determinantal representations for the 
correlation functions of the correspondent quantum model. 
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